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The  parabolic  equation  method  for  computing  acoustic  propagation  losses  was  used  to  model  the  results 
of  the  Tasman  Two  Sea  Tests.  In  these  tests  aircraft  dropped  SUS  (signals,  underwater  sound)  charges 
along  four  great  circular  paths  from  a hydrophone  near  South  Island.  New  Zealand,  in  westerly  directions 
across  the  South  Tasman  Sea  and  Southeast  Indian  Ocean.  The  data  were  processed  in  1/3  octave  bands 
from  16  to  1000  Hz  The  best  available  sound  speed,  bathymetric,  and  sediment  data  were  compiled  both 
from  Tasman  Two  Sea  Tests  data  and  from  historical  data.  These  data  were  used  to  make  parabolic  equa- 
tion runs  alone  the  radials  starting  at  the  stationary  receiver  for  63  and  1 25  Hz.  Great  care  was  required  in 
progressing  down  the  continental  slope  from  the  receiver  arid  over  sea  mounts  where  small  step  sizes  in  the 
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numerical  solution  were  required.  Accurate  modeling  of  sediments  was  also  required  in  these  areas.  Jn  the 
immediate  vicinity  of  the  hydrophone  a ray  theory  for  a wedge  gave  better  starting  values  for  the  parabolic 
method  than  did  a normal  mode  technique. 

Agreement  between  computations  and  experiment  was  generally  good  enough  to  substantiate  the 
techniques  used.  Poor  agreement  for  the  most  southerly  radial  suggested  inadequate  knowledge  of  the 
bathymetry  and  sediments  within  200  km  of  the  hydrophone.  Shadowing  behind  seamounts  was  found 
to  be  strongly  dependent  on  the  nature  of  the  sediment  cover. 


SUMMARY 

PROBLEM 

Apply  the  parabolic  equation  technique  to  model  the  acoustic  propagation 
characteristics  of  the  south  Tasman  Sea  and  southeast  Indian  Ocean  and  compare 
the  resulting  predictions  with  measured  acoustic  data.  Investigate  the  sensitivity 
of  the  predictions  to  vertical  and  horizontal  sample  step  sizes  and  to  the  specifica- 
tion of  the  initial  pressure  distribution  function.  Model  the  acoustic  shadow 
cast  by  a large  bathymetric  feature.  Study  the  acoustic  effects  of  horizontal 
changes  in  oceanography,  such  as  the  Subtropical  Convergence  and  the  Polar 
Front. 

RESULTS 

The  parabolic  equation  technique  was' used  to  help  interpret  Tasman  Two 
sea  test  data.  Computer  runs  were  made  with  realistic  sound-speed  and  bathy- 
metric data  for  ranges  up  to  2500  km  at  frequencies  of  63  and  1 25  Hz. 

A ray  trace  and  a normal  mode  prediction  were  used  to  initialize  the  para- 
bolic equation  model  in  a shallow,  coastal  environment,  hut  in  both  cases  the 
resulting  absolute  levels  predicted  in  deep  water  agreed  with  the  experimental 
levels  only  to  within  5 dB. 

Difficulties  were  encountered  in  attempting  to  determine  the  optimum 
vertical  sample  size  when  running  a parabolic  equation  prediction  from  shallow 
to  deep  water  and  vice  versa. 

The  gross  features  of  predicted  convergence  zones  (width  and  location)  as 
seen  from  the  coastal  site  agreed  well  with  the  experimental  data  out  to  ranges 
of  400  km. 

The  predicted  acoustic  shadow  of  substantial  underwater  features  was 
highly  sensitive  to  and  dependent  upon  the  assumed  properties  of  the  sediment 
cover.  The  observed  increase  in  shadow  depth  with  frequency  was  not  pre- 
dicted by  this  model. 

In  certain  range-dependent  environments,  the  signal  level  of  near-surface 
sources  increased  over  ranges  greater  than  1 000  km  and  the  “clarity”  of  the 
convergence  zone  structure  varied  dramatically  with  range. 
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RECOMMENDATIONS 


Study  the  use  of  the  parabolic  equation  method  in  a wedge-shaped  environ- 
ment to  establish  criteria  which  can  he  used  to  determine  the  optimum  values 
of  vertical  and  horizontal  sample  sizes  in  terms  of  water  depth,  frequency,  and 
bottom  slope,  with  the  aim  of  improving  the  accuracy  of  predicted  signal  levels. 

Apply  the  results  of  this  study  in  an  attempt  to  predict  the  correct  frequency 
dependence  of  the  acoustic  shadows  formed  by  large  underwater  features. 

Compare  the  efficiency  of  the  parabolic  equation  method  in  a wedge  with 
range-dependent  normal  mode  models. 

Continue  this  work  with  a quantitative  investigation  of  the  raw  predicted  data 
to  see  whether  the  variation  of  experimental  very  low  frequency  attenuation 
coefficients  can  be  partially  or  wholly  accounted  for  in  terms  of  gross  environ- 
mental characteristics  (sound-speed  profiles,  bathymetry,  sediment  types). 
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INTRODUCTION 

In  February  1975  the  Ocean  Science  Group  at  New  Zealand’s  Defence 
Scientific  Establishment  conducted  a very  low  frequency  acoustic  pro- 
pagation trial  (Project  Tasman  Two)  across  the  South  Tasman  Sea  into 
the  Southeast  Indian  Ocean  (Ref.  1 ).  The  trial  was  a joint  New  Zealand, 

United  States,  Australian  venture  involving  participation  of  additional 
personnel  from  the  Naval  Underwater  Systems  Center,  New  London,  and 
1 \N  Research  Laboratory,  Sydney.  The  aim  of  the  experiment  was 
ain  an  acoustic  “snapshot”  of  the  extensive  oceanic  region  to  the  south- 
New  Zealand  and  to  correlate  any  observed  changes  in  the  acoustic 
uion  characteristics  with  known  oceanographic  and  bathymetric 

lures. 

During  October  and  November  1975  the  authors  used  the  parabolic 
equation  technique  (Ref.  2)  to  make  extensive  propagation  loss  computa- 
tions along  the  experimental  paths.  The  purpose  of  these  computations 
was  to  determine  the  accuracy  with  which  the  experimental  results  could  be 
modeled  and  thus  determine  both  the  effectiveness  of  the  parabolic  equation 
technique  and  the  adequacy  of  the  environmental  data  in  describing  the 
propagation  conditions.  In  addition,  knowledge  could  be  gained  of  the 
propagation  mechanisms  revealed  by  the  data. 

The  environment  along  each  propagation  path  is  in  general  strongly  range 
dependent.  The  parabolic  equation  technique  is  known  to  work  well  in  a 
deep  ocean  situation,  where  the  only  significant  range-dependent  parameter 
is  the  sound-speed  profile.  However,  in  the  Tasman  Two  experiment  certain 
bathymetric  features  are  known  to  dominate  parts  of  the  environment.  Thus 
the  attempt  to  model  the  data  represents  a severe  test  of  the  parabolic  equation 
technique. 

The  experience  gained  in  this  use  of  the  parabolic  equation  technique  and 
the  increased  knowledge  of  its  capabilities  and  limitations  are  valuable  to 
NOSC’s  acoustic  prediction  abilities.  In  addition,  the  acoustic  properties  of 
the  ocean  area  included  in  these  experiments  and  computations  are  compara- 
tively unknown  and  these  results  are  a valuable  addition  to  the  data  base. 
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In  this  report  the  experimental  procedure  and  environmental  data  are  first  described. 
Next  the  parabolic  equation  technique  is  introduced,  with  considerations  of  step-sizes 
and  starting  values.  Finally,  the  particular  steps  of  this  investigation  are  outlined  and 
results  given.  Comparisons  between  various  computational  strategies  and  the  experi- 
mental results  are  shown. 


EXPERIMENTAL  DATA 
DESCRIPTION  OF  TEST  AREA 

The  hydrophone  system  used  in  Tasman  Two  was  bottom  mounted  at  a depth  of 
350  meters  near  45°S  167°E  (point  “T”)  and  cabled  back  about  6 miles  to  the  research 
ship  HMNZS  Tui  at  anchor  within  the  South  Island’s  Thompson  Sound.  RNZAF  air- 
craft flew  SUS  (signals,  underwater  sound)  drops  along  four  great  circle  paths,  or 
“radials,”  bearing  approximately  NW,  W,  WSW,  and  SW  from  the  hydrophone  site. 

The  SUS  were  set  to  detonate  at  1220  meters,  except  for  radial  WSW,  for  which  the  de- 
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F igure  L Approximate  location  of  principal  oceanographic  features  in  the  South  Tasman  Sea  and  Southeast 
Indian  Ocean.  Dashed  lines  indicate  acoustic  propagation  paths  investigated  during  Project  Tasman  Two. 
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tonation  depth  was  l)  14  meters.  Two  additional  drops  at  1 8 meters  were  made  along 
the  eastern  hall'  of  radials  NW  and  W.  The  aircraft  tracks  are  shown  in  f ig.  I . 

hollowing  the  SIJS  drops,  the  research  ship  I ni  made  an  oceanographic  cruise 
along  radial  VV  as  far  as  I55°l.  and  then  returned  to  New  Zealand  along  t he  complete 
length  of  radial  NW.  Bathymetry  was  recorded  continuously  and  velocity-tempera- 
ture-depth or  salinity-temperature-depth  lowerings  were  made  at  intervals  of  approx- 
imately 100  km  on  each  radial.  In  addition,  the  Australian  research  ship  HMAS 
Kimbla  undertook  an  oceanographic  cruise  over  the  large  bathymetric  feature  south 
of  Tasmania  known  as  the  South  Tasman  Rise.  The  oceanographic  and  bathymetric 
data  obtained  on  these  cruises  have  been  compiled,  together  with  relevant  historical 
data,  to  form  a fairly  complete  picture  of  the  environmental  conditions  along  each 
radial  (Ref.  3). 

Figure  I also  illustrates  the  main  oceanographic  features  of  the  region.  The 
oceanographic  conditions  in  the  area  have  been  summarized  by  (Jordon  ( Ref.  4), 
but  a brief  description  is  relevant  here:  The  Subtropical  Convergence,  which 
forms  the  boundary  between  the  relatively  warm,  saline  subtropical  water  to  the 
north  and  the  subantarctic  waters  of  the  south,  lies  in  a broad  band  approximately 
along  45°S.  From  an  acoustician’s  point  of  view  this  feature  marks  the  develop- 
ment ol  a relative  maximum  in  the  sound-speed  motile  at  a depth  ot  about  500 
meters.  This  maximum  is  fully  developed  at  the  southern  boundary  of  the  Sub- 
tropical Convergence  and  persists  to  the  south  until  it  disappears  at  the  Australasian 
Subantarctic  Front  near  50°S.  South  of  the  latter  feature,  the  sound-speed  mini- 
mum rises  steadily  from  an  initial  depth  of  about  1 250  meters  to  near-surtaee 
depths  at  about  55°S.  The  disappearance  of  the  deep  sound  channel  marks  the 
location  of  the  Polar  Front.  The  change  of  axis  depth  is  related  to  the  depth  ot 
the  salinity  minimum  produced  by  cold,  relatively  tresh  antarctic  intermediate 
water,  which  sinks  and  flows  northwards  at  the  Polar  front. 

ENVIRONMENTAL  DATA 

flic  oceanographic  and  bathymetric  conditions  along  the  three  radials  studied 
here  are  summarized  in  fig.  2.  3,  and  4,  as  tunctions  ol  distance  trom  point  I 
As  a general  rule,  the  continental  shelf  falls  away  steeply  from  the  New  Zealand 
coast  for  all  three  radials,  and  depths  of  4000  to  4500  meters  are  reached  within 
30  km. 

Radial  NW  (Fig.  2)  exhibits  the  most  uniform  bathymetry  once  clear  ol  the 
continental  slope,  flic  average  depth  is  about  4500  meters  and  there  are  no 
intrusions  above  3000  meters.  The  sound-speed  profiles  near  the  New  Zealand 
coast  are  typical  of  the  Subtropical  Convergence  in  summer,  with  surface  sound 
speeds  near  1510  m/s  and  an  inflection  point  at  about  500  meters.  I his  latter 
feature  is  present  only  in  profiles  measured  out  to  about  1000  km.  At  greater 
ranges  the  sound-speed  profiles  are  more  typical  of  warmer  subtropical  waters, 
as  indicated  by  the  rapid  increase  ol  surface  sound  speed,  which  reaches  a 
maximum  of  1 530  m/s  near  Sydney.  As  is  shown  later,  the  rather  subtle 
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Figure  2.  Radial  NW  showing  bathymetry  and  sound  speed  profiles.  The  tic  marks  indicate  the 
range  at  which  the  profile  was  measured  and  correspond  to  a sound  speed  of  1485  m/sec. 
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Figure  3.  Radial  WSW  showing  bathymetry  and  sound  speed  profiles.  The  tic  marks  indicate  the  range 
at  which  the  profile  was  measured  and  correspond  to  a sound  speed  of  1485  m/sec. 
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Figure  4.  Radial  SW  showing  bathymetry  and  sound  speed  profiles.  The  tic  marks  indicate  the  range 
at  which  the  profile  was  measured  and  correspond  to  a sound  speed  of  1485  m/sec. 


change  in  sound-speed  profile  seen  here  has  a marked  effect  on  the  predicted 
acoustics  of  the  region. 

The  bathymetry  of  radial  WSW  (Fig.  3)  is  dominated  by  the  South  Tasman 
Rise.  This  large  feature  is  approximately  500  km  across  at  the  base  and  rises 
to  a minimum  depth  of  850  meters.  Fxcept  in  the  vicinity  of  the  New  Zealand 
coast  the  oceanographic  conditions  along  this  radial  are  typical  of  Australasian 
subantaretic  waters.  The  surface  sound  speed  is  usually  less  than  1500  m/s 
and  the  profile  exhibits  a relative  maximum  near  500  meters. 

The  SW  radial  (Fig.  4)  has  several  interesting  features.  Four  different 
oceanographic  regions  can  be  identified  on  the  basis  of  sound-speed  profile. 

The  Subtropical  Convergence  region  is  seen  to  extend  out  to  about  450  km, 
after  which  subantaretic  conditions  prevail  until  the  Australasian  Subantaretic 
Front  is  reached  at  a range  of  approximately  1400  km.  The  channel  axis  then 
rises  over  the  following  200  km  until  the  deep  sound  channel  disappears  at 
about  1600  km.  Just  south  of  this  point  (the  Polar  Front)  the  surface  sound  speed 
drops  by  about  20  m/s  to  1465  m/s.  Apart  from  the  continental  slope,  the  bathy- 
metry is  dominated  by  two  major  features:  the  first  feature,  known  as  Resolution 
Ridge,  rises  to  a minimum  depth  of  2500  meters  at  a range  of  about  200  km.  The 
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second,  much  more  extensive  feature  forms  part  of  the  Southeast  Indian  Ridge 
(Ref.  5)  and  rises  about  2000  meters  at  a range  of  2000  km. 

ACOUSTIC  DATA 

The  received  shots  were  recorded  broadband  on  analog  tapes  which  were  returned 
to  the  Defence  Scientific  Establishment  for  analysis.  Each  signal  was  subsequently 
bandpass  filtered  through  1/3  octave  filters  at  1 6,  3 1 .5,  63, 125,  250,  500,  and 
1000  Hz.  After  the  signals  were  processed  by  an  analog  squarer-integrator  to 
obtain  the  energy  of  each  shot,  propagation  loss  versus  range  curves  were  produced 
for  each  frequency  band  using  source  levels  interpolated  from  curves  given  by 
Urick  (Ref.  6).  The  best  quality  data  were  obtained  in  the  125  anti  250  Hz  bands 
because  of  the  depth  of  the  shots.  The  two  lowest  frequency  bands  were  poor 
because  of  the  relatively  weak  source  level  and  high  ambient  noise  conditions 
at  this  particular  site.  The  unexpectedly  high  ambient  noise  level  was  thought 
to  be  associated  with  the  active  surf  conditions  experienced  near  the  hydrophones, 
which  were  placed  within  2 km  of  the  coastline. 

The  Tasman  Two  data  set  is  available  in  the  Naval  Sea  Systems  Command  ocean 
environmental  acoustic  data  bank  (NAVDAB)  system  (Ret.  7). 

SEDIMENTARY  STRUCTURES 

Detailed  data  on  the  properties  of  the  sediment  cover  in  this  area  are  sparse.  In 
general,  the  total  sediment  depth  is  found  to  be  of  the  order  100  meters  on  the 
continental  shell  and  South  Tasman  Rise  and  to  be  of  the  order  1000  meters  in 
the  ocean  basins.*  Several  boreholes  drilled  during  Leg  29  of  the  Deep  Sea  Drill- 
ing Project  (Ref.  X)  have  been  used  to  obtain  quantitative  spot  estimates  of 
sediment  composition.  In  addition,  I lout/,  compared  the  Deep  Sea  Drilling 
Project  data  with  seismic  reflection  profiler  results  and  thereby  extended  the 
coverage  away  from  the  immediate  vicinity  of  the  boreholes  (Ret.  9 and  10). 
following  Hamilton  ( Ref.  I I ) these  data  have  been  used  to  construct  geoacoustic 
models  of  the  sediment  structure  for  each  radial.  Details  of  the  models  are 
included  in  the  discussion  of  results. 

PROPAGATION  LOSS  CALCULATIONS 

SUMMARY  OE  PARABOLIC  EQUATION  METHOD 

I he  parabolic  equation  method  was  first  applied  to  underwater  acoustics 
problems  in  1973  by  Tappert  (Ref.  2).  although  a comprehensive  discussion 
of  the  technique  has  only  recently  been  published  ( Ref.  I 2).  The  essentials 
of  the  method  are  as  follows:  Sound  propagation  in  the  ocean  is  governed 


* Personal  communication  from  F.  J.  Davies,  Geophysics  Division,  197S. 
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hy  the  wave  equation  for  the  acoustic  pressure  excess,  p: 


&p  + k y n *(  r,z)p  = 0 (1) 

where  the  wavenumber  at  some  reference  sound  speed,  c , and  angular  frequency 
to,  is  given  by 

k()  - eo/c()  (2) 

and  the  refractive  index  is  given  hy 

n(r,z)  =cQ/c(r,z)  (3) 

Ihe  index  is  assumed  to  vary  only  with  depth,  r,  and  range,  r.  The  method  is, 
however,  applicable  to  eases  for  which  the  environment  varies  in  azimuth  as  well. 

II  n depends  on  depth  only,  Eq.  1 may  be  solved  as  a boundary  value  problem 
using  well-known  techniques  (Ref.  13).  The  solution  is  given  in  terms  of  a 
complete  set  of  cylindrical  waves ,pm(r,z),  known  as  “normal  modes.”  These 
waves  take  the  general  form 

Pm(r,z)  *Nm(kmr)/2  Um(z)e\p(ikmr)  (4a) 

where  Nm  is  a normalization  factor,  k is  the  horizontal  wave  number  and  U (z) 
gives  tlie  pressure  amplitude  as  a I unction  ol  depth.  I his  function  is  oscillatory 
and  has  m phase  reversals.  The  approximate  depth  increment  between  these  phase 
reversals  can  be  obtained  from  the  “vertical  wave  number,”  A'  . defined  by 

*m=^0-O'/2  <4ll> 

The  above  .approach  cannot  be  applied  to  range-dependent  environments. 
Instead,  one  of  several  approximate  methods  must  be  used.  The  simplest 
approximation  is  based  on  the  assumption  that  any  horizontal  gradients  are  too 
weak  to  induce  significant  mode  coupling;  the  modes  are  said  to  propagate 
ad  ia  bat  i ca  1 1 y ( Ret.  14).  This  means  that  acoustic  energy  initially  coupled  to  a 
given  mode  remains  carried  by  that  mode  as  it  adjusts  to  the  changing  environ- 
ment. This  method  may  be  implemented  by  recalculating  the  normal  modes  for 
each  new  environment  and  replacing  the  phase  by  its  integral  over  range.  This 
approach  is  ideally  suited  to  deep  ocean  situations  where  the  bathymetry  is 
unimportant.  However,  it  has  been  applied  in  modified  form  to  shallow-water 
situations  ( Ref.  I Sand  16). 

I he  basis  for  the  parabolic  equation  approximation  is  to  assume  that  the  total 
pressure  field,  p,  has  the  form 
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(5) 


-Zi 

p = U(r,z)r  e\p(ik0r) 

where  it  is  assumed  that  most  of  the  range  dependence  is  represented  hy  the 
second  term  and  that  U(r,z)  changes  only  slowly  with  range.  The  “parabolic” 
approximation  to  the  wave  equation  can  be  derived  by  substituting  Fq,  5 into 
Fq.  1.  A useful  form  of  this  approximation  is 

-2 ikn<L!l  dlu  + kl  \n(r,z)-\\U 
0 dr  dz2  ° 

which  is  valid  provided  the  following  conditions  hold: 
kQr>  1 

| EH  / 2 ikn<W  | <?  1 

3r2  dr 

The  first  condition  is  the  usual  for  field  approximation,  which  is  easily  satisfied 
in  the  ocean.  It  can  be  shown  that  the  second  condition  (Fq.  8)  amounts  to  a 
“small-angle”  approximation  (Appendix).  Fitzgerald  ( Kef.  17)  has  shown  that 
Fq.  6 is  valid  for  a small  cone  of  angles  which  may  he  oriented  in  any  direction. 
The  mean  angle  is  determined  by  c0,  and,  in  order  to  minimize  (he  error,  it  is 
usual  to  select  a value  midway  between  the  maximum  and  minimum  sound  speeds 
encountered. 

These  approximations  convert  the  boundry  value  problem  to  an  initial  value 
problem.  In  other  words,  given  the  function  U at  some  range,  rQ , Fq.  6 can  be 
used  to  “march”  the  solution  out  to  ranges  greater  than  r(J.  This  technique  can 
be  implemented  by  making  use  of  the  so-called  “split-step  Fourier  algorithm” 
developed  by  Tappcrt  and  Hardin.  A detailed  derivation  is  given  in  Ref.  12  and 
we  quote  only  the  result  here: 

U(r  + Ar.z)  ^exp  [ikQ(n2-  1)1  F’  ' ) Fs[  (/(r,z)j^  <’h 

The  symbols  Fs  and  F'1  denote  the  Fourier  transform  into  “S”  space  and  its 
inverse,  and  A r is  the  range  increment. 

The  following  comments  hold  regarding  the  usefulness  of  l.q.  l>  anti  the 
parabolic  equation  method  in  general: 

Accuracy,  liquation  l>  is  an  exact  solution  to  Fq.  6 (which  itself  is  an  approxi- 
mation to  the  wave  equation)  only  if  the  refractive  index  is  a constant  ( Kef.  I 2). 
However,  vertical  changes  of  sound  speed  are  typically  less  than  57  and  strong 
near-surface  horizontal  sound-speed  gradients  are  of  the  order  10  1 m/s/km 
(Kef.  14).  On  this  basis  it  is  assumed  that  Fq.  9 gives  a valid  solution  to  the 
wave  equation  for  a real  ocean. 


(6) 

(7) 

(8) 
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Density.  The  parabolic  equation  method  does  not  take  the  density  change 
at  the  water-sediment  interlace  into  consideration.  Thus,  a somewhat  greater 
bottom  reflection  loss  than  normal  is  to  be  expected.  However,  numerical 
results  indicate  that  the  difference  is  small.  Tor  example,  a density  contrast 
of  1.25: 1 causes  a decrease  in  reflection  loss  of  less  than  0.5  d B/bounce  at  low 
grazing  angles.  In  contrast,  a 5"  jump  in  sound  speed  at  the  sediment  inter- 
face causes  a 4 to  5 dB  decrease  for  the  same  case  (Tig.  5).  Thus,  the  inability 
to  include  density  in  the  model  is  unlikely  to  be  a source  of  serious  error. 

Sampling  Rates.  The  numerical  solution  of  Tq.  9 relies  on  the  use  of  the 
fast  T'oiifiei  transform  to  advance  the  solution  in  range  at  a rapid  rate.  However 
to  avoid  error  care  must  be  taken  to  sample  the  pressure  field  with  a vertical 
step-size,  Az,  which  is  ideally  about  1/5  the  minimum  vertical  “wavelength” 
present.  This  sets  a practical  upper  limit  to  the  number  of  modes  (and  thus  the 
frequency)  that  can  be  handled,  because  the  computer  run-time  increases  as  n 
\og(n)  where  n is  the  number  of  fast  Fourier  transform  points.  Similarly  the 
horizontal  step-size,  Sr,  is  limited  by  the  frequency  and  magnitude  of  horizontal 
gradients.  However  there  appears  to  be  no  deterministic  criteria  which  can  be 


Figure  5.  Bottom  reflection  loss  obtained  from  normal  mode  program  and  converted  to  units  of  dB/bounce. 
Water  depth  350m.  Frequency  63Hz. 
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applied  to  establish  Ar.  Usually  a trial  and  error  approach  is  taken  in  which 
Ar  is  progressively  reduced  until  the  solution  converges.  This  implies  that  in 
the  interests  of  efficiency  parabolic  equation  computer  models  need  to  optimize 
dynamically  both  vertical  and  horizontal  step  sizes  according  to  the  local  environ- 
mental conditions. 

Phase  Error,  liquation  6 can  be  used  to  illustrate  the  propagation  of  a single 
mode  under  the  parabolic  equation  approximation  (Ref.  18).  It  is  found  that 
the  mode  propagates  as: 

p,n*HmUm  (z)  (kmr)Vl  exp[/(fr0-**  / 2k0)r)  (10) 

Comparing  l.q.  10  with  l.q.  4a  and  4b,  it  is  apparent  that  whereas  the  mode 
amplitude  is  unaffected,  the  parabolic  equation  phase  is  a first-order  approxi- 
mation of  the  true  phase,  (Ay  - )'Ar.  t hus  in  the  parabolic  equation  approxi- 

mation phase  sensitive  effects,  such  as  l.loyd  Mirror  interference  patterns,  develop 
a range-dependent  offset  relative  to  a normal  mode  prediction. 

Reciprocity.  It  is  well  known  ( Ref.  1 3 ) that  the  normal  mode  solution,  which 
is  valid  for  environments  with  no  horizontal  gradients,  is  invariant  under  the 
exchange  of  source  and  receiver  coordinates.  This  phenomenon  is  known  as 
“reciprocity.”  When  the  medium  is  range  dependent  it  is  still  theoretically  true 
that  reciprocity  holds.  Experience  with  ocean  fronts  indicates  that  reciprocity 
does  hold  to  a good  approximation  in  deep  water,  even  in  regions  ot  extreme 
horizontal  gradient.*  I 'here  fore,  reciprocity  is  assumed  to  hold  lor  the  radials  in 
in  Tasman  Two.  Accordingly  all  predictions  are  for  a bottomed  source  radiating 
from  the  hydrophone  site  to  a towed  receiver  at  the  shot  depth.  This  greatly 
reduces  the  computer  runtime  because  the  only  alternative  is  to  calculate  the 
propagation  loss  back  to  the  hydrophone  from  each  point  of  interest.  In  fact, 
some  differences  between  up-slope  and  down-slope  calculations  have  been  observed 
and  are  discussed  in  the  next  section. 

Cost.  | he  time  required  to  make  parabolic  equation  runs  is  known  to  be 
high  and  can  be  a limiting  factor  in  its  use.  The  longer  runs  reported  in  this 
report  required  roughly  I / 2 to  1 hour  each  on  the  Univae  I 1 10  at  a cost  on 
low  priority  time  of  $200  to  $400.  The  longest  times  were  required  for  those 
runs  over  the  South  Tasman  Rise,  which  required  small  horizontal  step  size,  as 
is  discussed  elsewhere. 

PROGRAM  CONSIDERATIONS 

The  parabolic  equation  model  used  was  a modified  version  of  an  existing 
NOSC  computer  program.  The  modified  program  accepts  an  unlimited  number 
of  input  sound-speed  profiles  at  specified  ranges  and  linearly  interpolates  the 
refractive  index  between  the  two  profiles  which  straddle  the  current  range. 
Similarly,  the  bathymetry  is  represented  as  a sequence  ol  linear  segments 
between  the  input  depth-range  data  points. 

* Personal  communication  from  M.  A.  Pedersen  and  D.  I .Gordon,  1975. 
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F I lie  magnitude  of  both  the  range  increment,  Ar,  and  the  jump  discontinuity 

in  sound  speed  at  the  sediment  interlace  can  also  be  specified  as  a function  of 
range.  However  the  specification  of  the  bottom  model  (sediment  type  and 
thickness)  is  a constant  for  each  run.  If  necessary,  the  bottom  model  can  be 
changed  by  first  generating  a new  initial  pressure  distribution  at  a given  range 
and  then  rerunning  the  program  to  greater  ranges  with  the  new  sediment 
structure.  Similarly,  the  depth  increment  A z cannot  be  changed  without  re- 
starting the  program. 

Initial  pressure  distribution  functions,  U(r^,z)  (which  are  referred  to  as  “starters”), 
were  calculated  using  both  normal  modes  anti  rays.  In  deep  water  it  is  usual  to 
start  the  parabolic  equation  model  with  a (iaussian-shaped  initial  pressure  distri- 
bution close  to  the  source.  1 he  Gaussian  shape  is  used  because  it  filters  out  high- 
angle  energy  t. which  would  be  lost  at  short  range  anyway)  and  prevents  aliasing 
of  this  energy  into  low  angles  during  the  fast  Fourier  transform  sampling  process. 
However,  this  technique  was  not  suitable  to  start  the  Tasman  Two  predictions 
because,  under  the  assumed  conditions  of  reciprocity,  the  “source”  was  placed  on 
the  bottom  in  shallow  water  at  the  top  of  the  continental  shelf. 

The  mode-generated  starters  involve  two  approximations  (Tig  6a):  ( I ) it  is 
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necessary  to  assume  the  water  depth  is  constant  at  the  top  of  the  continental 
slope  in  order  to  apply  mode  theory;  and  (2)  numerical  difficulties  impose  an 
upper  limit  of  about  60  degrees  to  the  equivalent  ray  angle.*  The  latter  approxi- 
mation restricts  the  minimum  range  at  which  the  mode  starter  is  accurate  to  about 
200  meters,  t hus,  the  bathymetry  was  held  at  a constant  depth  of  350  meters 
tor  a distance  200  meters  down  range  from  the  source. 

The  ray  starter  was  used-in  an  attempt  to  reduce  these  approximations.**  A 
constant  bottom  slope  appropriate  to  the  radial  concerned  was  applied  near  the 
source  (l  ig  6b),  and  the  pressure  distribution  at  I km  was  calculated  assuming 
a constant  1.5  dB/hounce  rellection  loss.  'This  figure  was  based  on  the  results 
of  mode  calculations  for  various  bottom  models.  Rays  in  the  forward  direction 
up  to  ± 90  degrees  where  included. 

It  is  ol  course,  necessary  to  consider  very  high  angle  energy  close  to  the  source, 
because  the  angle  is  considerably  reduced  during  propagation  down  the  slope. 
Ignoring  the  intermode  coupling  caused  by  the  bathymetric  gradient,  the  minimum 
requirement  is  that  energy  which  insonilies  conveyance  /.ones  when  it  reaches 
d«>ep  water  must  be  treated  properly  at  the  top  of  the  slope.  On  this  basis,  the 
vertical  sample  interval  may  be  expressed  as 

A z ~ 2/1 /Sin  (||) 

where  //  is  the  water  depth  and  in  denotes  the  highest  order  mode  required  in 
deep  water,  f or  example,  at  1 00  11/  it  is  necessary  to  consider  about  1 50  modes 
in  deep  water  to  include  some  bottom  bounce  energy,  l or  a water  depth  of 
4500  meters,  f.q.  I I indicates  that  a step  size  ol  I 2 meters  is  adequate.  However 
at  the  lop  ot  the  continental  slope  where  //  is  only  350  meters,  the  step  size  is 
reduced  to  about  I meter. 


DISCUSSION 


PROCEDURES  AND  FINDINGS 

In  view  of  the  anticipated  difficulties  associated  with  initializing  the  parabolic 
equation  model  on  the  continental  slope,  it  was  decided  to  attempt  predictions 
ol  the  average  signal  level  over  the  radials  before  investigating  specific  details 
ot  the  acoustic  data.  Consequently  the  initial  effort  was  directed  towards  a 
study  ol  the  effect  on  absolute  level  of: 

1.  The  vertical  and  horizontal  step  sizes 

2.  The  sedimentary  structure  of  the  bottom 

3.  The  specification  of  starter  functions 

Subsequent  work  reported  here  discusses  the  convergence  zones  of  radial  NW.  the 
acoustic  shadow  cf  the  South  Tasman  Rise,  and  some  distinctive  acoustic  features 
related  to  changes  in  the  sound-speed  profile. 

* Personal  communication  from  D.F.  Gordon.  1975. 

**  Personal  communication  from  R.  White  ami  D.F. Gordon,  1975. 


GEOACOUSTIC  MODELS 

Various  geoacoustic  models  were  investigated  before  selecting  the  one  we  refer 
to  as  a sediment  bottom  (Fig.  7).  The  constraints  which  led  to  this  model  being 
adopted  were: 

1.  In  the  absence  of  firm  data  on  the  sedimentary  structure  of  the  continental 
shell,  the  model  should  not  be  too  detailed. 

2.  The  model  should  utilize  accepted  parameters  to  specify  the  acoustic  pro- 
perties of  sediments,  such  as  those  given  by  Hamilton  (Ref.  I I ). 

3.  In  view  of  the  inability  to  change  the  bottom  model  with  range,  the  model 
should  be  capable  of  describing  the  deep  ocean  basins. 

4.  The  model  should  be  compatible  with  the  requirements  of  the  normal  mode 
program  used  to  generate  starters. 
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The  bottom  was  therefore  represented  as  a series  of  fluid  layers,  in  each  of  which 
the  sound  speed  varies  as 

c(  z)=ci  [ 1 -27(.(  z-z.  ) / ct  J-'I 2 

where  C.  and  7.  are  the  sound  speed  and  sound-speed  gradient  at  the  top  of  the 
/th  layer  and  z . is  the  layer  depth,  liquation  12  is  useful  in  that  the  normal  modes 
can  he  expressed  in  terms  ol  Airy  functions,  hut  is  limited  in  that  the  derived 
sound  speed  and  gradient  take  on  physical  values  over  a restricted  depth  range 
only.  Accordingly  was  set  to  the  measured  value  at  point  “T”  and  the  remain- 
ing parameters  specified  as  in  I able  1 . 1 he  layered  structure  has  no  direct  experi- 
mental significance,  being  imposed  to  ensure  a smooth  progression  of  sound  speed 
and  gradient  with  depth,  within  the  limitations  of  l.q.  I 2,  and  in  addition,  to 
ensure  that  solutions  for  many  high-order  modes  could  be  found.  The  negative 
gradient  in  the  5th  layer  is  necessary  to  satisfy  the  mechanics  of  the  mode  program 
and  has  no  effect  on  the  solution,  as  it  is  at  a depth  well  below  that  to  which  acoustic 
energy  penetrates  at  a significant  intensity. 


Table  1.  Bottom  layer  parameters  for  the  sediment  model. 


I o compensate  for  the  inability  to  include  density  in  the  model  and  to  allow 
for  greater  bottom  reflectivity  on  the  continental  slope  as  compared  with  the 
ocean  basins,  a variable  jump  discontinuity  was  applied  at  the  water-sediment 
interface.  Usually  this  jump  was  taken  as  11'-'  on  the  slope  but  was  tapered  to 
zero  on  reaching  the  ocean  basins. 

f inally,  the  attenuation  coefficient  in  the  bottom,  a,  was  taken  as  a constant 
0. 1 l'uz  dB/km,  corresponding  to  a porosity  of  65'a  or  greater  (sand-silt-clay  or 
liner  grain  sediment ).  However  because  of  the  jump  discontinuity  allowed  at 
the  surface,  it  is  recognized  that  this  approach  is  somewhat  unphysical  in  view 
of  the  well-known  empirical  relationships  between  sound  speed,  porosity,  and 
attenuation  ( Ref.  I 1 ). 
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ABSOLUTE  LEVELS 

To  begin  the  study  of  absolute  signal  levels,  we  investigated  the  use  of  the 
parabolic  equation  model  in  an  environment  similar  to  that  near  the  hydrophone 
site.  Eli  is  was  modeled  as  a sound  channel  having  a constant  depth  of  350  meters 
overlying  a sediment  bottom.  The  sound-speed  profile  was  as  measured  at  point 
“I  ”,  f igure  8 presents  a comparison  of  the  normal  mode  prediction  and  the 
parabolic  equation  prediction.  The  latter  was  initialized  by  the  mode  solution 
at  200  meters  and  marched  out  with  a range  increment  of  10  meters  and  a 
depth  increment  of  I 2 meters.  Apart  from  the  phase  difference  of  about  500 
meters,  the  predictions  are  very  close  at  ranges  greater  than  about  2 km.  Investi- 
gation of  the  mode  solution  showed  that  fewer  than  10  modes  contributed  at 
these  ranges,  in  which  case  the  12-meter  sample  interval  is  expected  to  give  an 
accurate  solution.  At  short  ranges  where  many  modes  contribute,  the  agree- 
ment is  not  as  good. 

fuming  to  the  experimental  situation,  we  consider  radial  NW,  which  has 
the  simplest  environmental  conditions  at  short  range.  Figure  0 summarizes 
tlie  best  results  obtained  using  the  sediment  bottom  model  at  63  Hz.  Three 


Figure  8.  Comparison  of  predicted  propagation  loss  in  shallow  water  channel  overlying  sediment  bottom. 
Water  depth  350  meters,  frequency  63  Hz.  No  density  change  at  sediment  interface. 
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cases  are  illustrated:  curve  a was  derived  usin  ',  a normal  mode  starter  at  200 
meters  with  A z initially  6 meters  hut  doubled  to  I 2 meters  at  a range  of  20 
km  (water  depth  2000  meters).  Ar  was  10  meters  down  the  shelf  and  200 
meters  in  deep  water.  An  I 1%  jump  in  sound  speed  at  the  sediment  interlace 
was  applied  on  the  shelf  out  to  a range  of  27  km.  This  was  tapered  to  zero 
over  the  following  3 km.  Propagation  loss  was  calculated  at  2-km  intervals 
hut  is  presented  as  a 5-point  average  in  order  to  smooth  the  curve. 

Comparison  of  this  curve  with  the  experimental  data  shows  that  the  pre- 
dicted level  is  about  5 dB  lower  than  the  measured  I vel  Although  this  could 
be  due  to  a source-level  error,  we  observe  that  the  prediction  obtained  using 
a 1 2-meter  vertical  sample  throughout  (curve  hi  is  significantly  different  in 
detail.  This  is  an  indication  that  the  solution  has  not  stabilized,  most  probably 
because  Az  needs  to  be  reduced  further  Assuming  a'  •>  t <■()  significant  modes 
in  deep  water  at  63  Hz,  liq.  I 1 indicates  that  A.  si  . os  than  3 meters 

at  the  top  of  the  shelf.  Unfortunately  it  was  m i p s mn  the  existing  model 

for  Az  this  small  without  further  modification  to  tu  ; m 


Figure  9.  Radial  NW.  Comparison  of  parabolic  equation  propagation  loss  predictions  with  experimental 
data  at  63  Hz. 
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However  checks  were  made  on  some  of  the  other  parameters  of  interest:  To 
investigate  the  effect  of  different  sediment  structure,  the  prediction  was  run  down 
the  slope  with  a bottom  modelled  by  a sediment  layer  100-meters  thick,  porosity 
65'/,  attenuation  coefficient  OI*Hz  and  sound  speed  jump  3.3%,  overlying  base- 
ment rock.  The  predicted  level  (not  shown)  was  about  5 dB  lower  than  curve  a, 
presumably  because  of  increased  absorption  within  the  sediment  due  to  the  reduced 
bottom  reflectivity.  It  is  interesting  that  the  presence  of  the  strongly  reflecting 
layers  of  basement  rock  (sound  speed  5000  m/s)  did  not  raise  the  signal  level. 

This  indicates  that  little  acoustic  energy  penetrates  as  far  as  the  basement,  which 
is  at  a depth  of  4 wavelengths  at  this  frequency,  and  to  some  extent  justifies  Ihe 
use  of  the  sediment  model  on  the  continental  shelf. 

In  addition,  a prediction  was  made  with  the  ray  starter  which,  as  discussed 
earlier,  uses  a constant-slope  approximation  to  the  bathymetry  near  “T”  and 
includes  all  initially  forward-going  energy  (curve  c).  The  sampling  intervals  and 
all  other  parameters  are  the  same  as  for  curve  a.  We  observe  that  because  of  the 
greater  water  depth  at  I km  (600  meters),  where  the  parabolic  equation  starter 
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Kgure  10.  Radial  NW:  Comparison  of  parabolic  equation  propagation  loss  prediction  with  experimental 
data  at  63Hz. 
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is  generated,  the  6-meter  vertical  sample  interval  is  a better  approximation  in 
this  case.  Comparing  curve  c with  curve  a,  the  ray  starter  yields  a predicted 
level  about  9 dB  greater  than  the  mode  starter.  However,  the  overall  level 
exceeds  that  of  the  experimental  data  by  a few  dB.  A possible  cause  of  this 
excess  level  is  the  simple  1.5  dB/bounce,  constant  reflection  loss  model  assumed 
for  the  ray  calculation.  I rom  fig.  5 this  low  loss  would  tend  to  exaggerate  the 
high-angle  energy  content. 

figures  10  and  1 I illustrate  the  predictions  for  deep  shots,  over  the  complete 
span  of  radial  NW,  at  63  11/  and  I 25  Hz,  respectively.  The  parabolic  equation  model 
was  initialized  with  the  ray  starter  in  both  cases  and  apart  from  frequency,  all 
parameters  are  as  lor  curve  c in  fig.  9.  These  figures  represent  the  best  fit  to  the 
overall  signal  levels  on  this  radial.  However,  in  view  of  the  previous  discussion  it 
is  expected  that  the  degree  of  agreement  would  be  significantly  improved  by  the 
use  of  smaller  Ac  over  the  first  few  kilometers  of  range  and  a better  reflection 
loss  model  when  calculating  the  starter. 


o 


figure  11.  Radial  NW:  Comparison  of  parabolic  equation  propagation  loss  prediction  with  experimental 
data  at  1 25  Hz. 
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CONVERGENCE  ZONES 

The  acoustic  data  for  the  shallow  shots  on  radial  NW  are  interesting  in  that 
the  convergence  zone  structure  is  more  pronounced,  both  in  the  higher  frequency 
bands  and  at  longer  ranges  (Ref.  I ).  Figures  I 2 and  1 3 illustrate  the  convergence 
zone  predictions  ( 18-meter  shots)  for  this  radial  at  63  and  1 25  Hz,  respectively. 
These  predictions  were  generated  at  the  same  time  as  those  in  Fig.  10  and  1 I . 

The  bathymetry  has  been  superimposed  on  these  figures  to  indicate  the  scale  of 
the  continental  slope,  which  drops  down  to  the  deep  ocean  within  the  first 
convergence  zone. 

Comparing  the  experimental  and  predicted  data  at  63  Hz,  we  see  that  the 
experimental  zones  are  best  defined  at  ranges  greater  than  300  km  but  are  about 
half  a zone  ahead  of  the  prediction  in  phase.  At  125  Hz  ( Fig.  13),  the  measured 
zones  over  the  first  300  km  arc  better  defined  than  at  63  Hz  and  compare  quite 
well  with  the  prediction.  However  as  in  Fig.  1 1 , the  predicted  level  is  too  high 
on  the  average.  At  longer  ranges  the  phase  difference  seen  in  Fig.  1 2 is  evident  here 
also. 
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Figure  12.  Radial  NW:  Comparison  of  predicted  convergence  zone  propagation  loss  with  experimental 
data  at  63  Hz. 
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Figure  13.  Radial  NW:  Comparison  of  predicted  convergence  zone  propagation  loss  with  experimental 
data  at  1 25  Hz. 


The  most  significant  difference  between  the  data  and  the  predictions  is  that 
the  predicted  level  drop  between  peak  and  shadow  is  about  twice  that  measured. 
Acoustic  energy  is  of  course  carried  into  the  shadow  /.one  by  bottom-reflected 
modes.  Therefore,  the  lack  of  agreement  may  point  to  a need  for  a bottom 
model  with  lower  reflection  loss.  It  is  also  possible  that  sufficient  bottom- 
reflecting  modes  were  not  included  in  the  calculation  because  the  vertical 
sample-size,  A z,  was  too  coarse.  A third  possible  reason  for  the  discrepancy 
is  that  we  are  comparing  a continuous  wave  prediction  with  data  derived  from 
a I /3-octave  energy  analysis.  However,  any  smoothing  induced  by  the  finite- 
band  analysis  is  minimal  in  this  case,  because  the  gross  features  of  convergence 
zones  are  essentially  frequency  independent.  In  view  of  the  operational  signifi- 
cance of  these  effects  in  areas  where  refracted-surfacc  reflected  ( RSR)  propaga- 
tion is  marginal,  it  was  unfortunate  that  the  investigation  could  not  be  pursued 
further. 

ACOUSTIC  SHADOWING 

On  radial  WSW  the  SIJS  were  set  to  detonate  at  a depth  of ‘)|4  meters,  which 
is  just  above  the  axis  of  minimum  velocity  but  slightly  below  the  peak  of  the 
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figure  14.  Radial  WSW:  Experimental  data  showing  acoustic  shadow  of  South  Tasman  Rise.  Reproduced 
with  permission  from  Ref.  L 


South  I asman  Rise.  I bis  shoals  to  a minimum  depth  of  850  meters  (Fig.  3). 

I lie  principal  feature  of  the  acoustic  data  for  this  radial  is  the  “shadow”  cast 
hy  the  Rise,  which  of  course  does  not  act  as  a complete  stop,  because  some  acoustic 
energy  is  channeled  over  the  top  via  bottom-reflected  paths.  By  extrapolating  the 
propagation  lose  data  east  of  the  Rise  back  into  the  shadow  region,  the  depth  o'" 
the  acoustic  shadow  can  be  measured.  The  shadow  at  63  11/  is  close  to  7 dB  (above 
cylindrical  spreading  loss)  and  it  increases  at  a rate  of  about  4.5  dB 'octave  over 
the  higher  frequency  bands  (I'ig.  14).  The  data  in  the  16- and  31.5-11/  bands  are 
of  insufficient  quality  to  obtain  accurate  estimates. 

In  order  to  reduce  the  computer  run-time  for  this  radial,  two  approximations 
were  made,  f irst,  the  parabolic  equation  program  was  initiated  at  l)0()  km  using 
a normal  mode  starter.  Consequently  the  absolute  levels  predicted  in  this  section 
arc  somewhat  arbitrary.  Second,  a fixed  sound-speed  profile  (measured  at  1470 
km  in  l ig.  3)  was  used  on  the  assumption  that  the  bathymetry  would  dominate 
the  shadowing  process.  Therefore  the  only  environmental  variables  here  are  the 
bathymetry  and  bottom  model,  l or  each  run  the  horizontal  step-size,  Ar,  was 
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Figure  IS.  Radial  WSW:  Comparison  of  parabolic  equation  propagation  loss  prediction  with  data  at  63  Hz. 

Sediment  bottom  model  with  no  sound-speed  jump  at  sediment  interface. 

progressively  reduced  ironi  200  meters  in  deep  water  (at  63  Hz),  to  20  meters  at 
the  top  of  the  Rise.  At  1 25  Hz  Ar  was  half  these  values.  However,  the  vertical 
sample-size,  A z,  was  fixed  at  1 2 meters  for  all  cases. 

To  illustrate  the  importance  of  the  bottom  model  in  this  situation  Tig.  1 5 gives 
the  prediction  obtained  at  63  llz  using  the  sediment  bottom  <1  ig.  7)  with  no  sound 
speed  jump  at  the  water-sediment  interface.  It  is  apparent  that  the  feature  acts  as 
an  almost  complete  acoustic  stop,  most  of  the  energy  being  absorbed  within  the 
sediments.  The  bottom  loss  was  reduced  considerably  by  imposing  a jump  in  sound 
speed  at  the  sediment  interface.  The  magnitude  of  the  jump  was  tapered  from  zero 
either  side  of  the  Rise  to  I \7<  at  depths  less  than  1500  meters.  At  63  Hz  this 
bottom  model  gives  a good  fit  to  the  data  (solid  line  in  lug.  16),  but  it  is  interest- 
ing that  at  125  llz  (I  ig.  1 7)  the  shadow  obtained  using  the  same  model  is  not 
sufficiently  deep.  The  predicted  drop  in  level  west  of  the  South  Tasman  Rise 
at  125  11/  is,  in  fact,  5 dlt  less  than  that  at  63  llz.  The  gap  in  predicted  level 
was  caused  by  the  source  (or  receiver)  depth  intersecting  the  sediment  at  the 
top  of  the  sea  mount. 
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Figure  16.  Radial  WSW:  Comparison  using  sediment  models  with  sound-speed  jumps;  frequency  63  Hz. 


In  an  attempt  to  improve  the  prediction,  we  used  a more  detailed  bottom  model 
of  tlie  Rise  based  on  seismic  surveys  of  the  area  and  data  logged  from  borehole  281. 
DSDP  Leg  29,  which  was  drilled  on  the  south  flank  of  the  South  Tasman  Rise 
(Ret.  8.  9,  and  10).  I he  model  (t  ig.  1 8)  consists  of  a 25-meter  layer  of  strongly 
absorptive,  very  fine  sand  overlying  I 75  meters  of  less  absorptive,  sand-silt-clay. 

(I  or  this  reason  it  is  referred  to  as  the  “two-layer”  model.)  Basement  occurs  at 
200  meters  in  the  model  and  imposes  a 300  - jump  in  sound  speed  (Table  2). 


Table  2.  Bottom  layer  parameters  for  the  two-layer  model. 
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Figure  17.  Radial  WSW:  Comparison  using  sediment  models  with  sound-speed  jumps,  frequency  1 25  Hz. 

To  implement  the  necessary  changes  to  the  bottom  model,  the  program  was  run 
in  three  stages.  First  the  prediction  was  run,  using  the  sediment  model,  from  900 
to  1280km,  at  which  point  a new  starter  was  stored.  The  prediction  was  then 
run  over  the  South  Tasman  Rise,  using  the  two  layer  model,  out  to  I860  km,  where 
a second  starter  function  was  stored.  Finally  the  prediction  was  continued  using 
the  sediment  model  again. 

I'he  results  are  illustrated  by  the  broken  lines  in  Fig.  16  and  17.  We  sec  that 
the  predicted  shadow  at  1 25  Hz  has  deepened  and  more  closely  fits  the  data. 
However,  the  shadow  at  63  II/.  is  now  considerably  deeper  than  the  observed 
shadow. 

These  results  arc  encouraging  in  that  both  the  onset  of  the  shadow  and  the 
initial  rate  of  decrease  in  signal  level  behind  the  Rise  agree  well  with  the  data. 

At  the  same  time,  the  results  are  disappointing  in  that  the  predicted  frequency 
dependence  of  the  acoustic  shadow  is  opposite  to  that  observed:  l or  both  bottom 
models  the  level  drop  behind  the  Rise  decreases  with  increasing  frequency.  This 
discrepancy  could  be  related  to  the  choice  of  vertical  sample  size,  Ar,  which 
remained  constant  at  I 2 meters  throughout.  At  the  top  of  the  Rise,  where  the 
water  depth  is  850  meters,  a sample  size  of  about  3 meters  would  be  more 
appropriate  at  125  Hz,  assuming  120  modes  are  important  (Fq.  1 1). 
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Figure  18.  Important  parameters  of  “two-layer”  bottom  model. 

OCEANOGRAPHIC  EFFECTS 

The  acoustic  propagation  data  from  radial  SW  is  dominated  by  the  presence  of  a 
high-loss  region  at  short  ranges.  After  the  effects  of  cylindrical  spreading  are  removed, 
the  signal  levels  at  ranges  greater  than  450  km  are  about  14  dB  lower  than  the 
levels  extrapolated  from  about  150  km  (Fig.  Ilf).  In  contrast  to  the  frequency- 
dependent  level  drop  seen  west  of  the  South  Tasman  Rise,  the  level  drop  on  radial 
SW  appears  to  be  independent  of  frequency.  South  of  this  anomaly,  two  distinct 
attenuation  regimes  can  he  defined.  The  first  regime  extends  from  the  acoustic 
anomaly  out  to  a range  of  about  1 200  km  and  apparently  coincides  with  the  sub- 
antarctic  water  mass  lying  between  the  Subtropical  Convergence  and  the 
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Figure  19.  Radial  SW:  Experimental  data  showing  acoustic  anomaly  for  ranges  <450  km  and  high  attenuation 
zone  south  of  Polar  Front.  Reproduced  with  permission  from  Ref.  L 


Australasian  Subantarctic  Front.  The  second  regime  extends  southwest  from 
the  Australasian  Subantarctic  Front  over  the  remainder  of  the  radial. 

Figure  20  illustrates  the  parabolic  equation  prediction  at  125  Hz  for  this  radial 
using  the  ray  starter.  The  sediment  bottom  model  was  used  over  the  entire  radial, 
with  an  I 1%  jump  in  sound  speed  at  the  sediment  interface  on  the  continental 
slope.  The  vertical  sample-size  was  initially  6 meters,  but  was  increased  to  12 
meters  for  ranges  greater  than  20  km.  In-shore  of  the  Resolution  Ridge  at  ranges 
less  than  200  km  the  prediction  follows  the  trend  of  the  data  but,  as  discussed 
for  radial  NW,  the  average  level  is  a few  decibels  too  high.  However,  just  south 
of  the  Ridge,  where  the  experimental  level  is  anomalously  low.  the  prediction 
shows  no  evidence  of  a corresponding  level  drop.  The  reasons  for  Ibis  are  not 
clear.  One  possibility  is  that  the  Resolution  Ridge  rises  to  much  shallower  depths 
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Figure  20.  Radial  SW:  Comparison  of  propagation  loss  predictions  with  experimental  data  at  1 25  Hz. 


than  indicated  by  bathymetric  charts  of  the  area.  However,  the  lack  of  a fre- 
quency dependence  similar  to  the  shadow  of  the  South  Tasman  Rise,  raises  the 
interesting  possibility  that  this  anomaly  is  related  to  some  oceanographic  pro- 
perty of  the  area  rather  than  to  the  bathymetry.  The  fact  that  there  is  no  notice- 
able level  drop  associated  with  the  Southeast  Indian  Ridge,  a massive  feature 
some  350  meters  higher  than  the  Resolution  Ridge,  lends  support  to  this  hypo- 
thesis. In  any  case,  it  is  clear  that  the  environment  near  the  Resolution  Ridge 
has  not  been  adequately  represented. 

A second  interesting  feature  of  this  prediction  is  the  change  in  character  that 
occurs  near  1 500  km.  Evidently  this  is  caused  by  the  disappearance  of  the  deep 
sound  channel  at  the  Polar  Front  (Fig.  4).  However,  south  of  the  Polar  Front 
the  prediction  does  not  show  a marked  increase  in  slope  as  docs  the  experimental 
data.  It  is  possible  that  this  increase  in  attenuation  is  caused  by  rough  surface 
effects  not  included  in  the  model:  Presumably  any  surface-induced  coupling  of 
energy  into  bottom  bounce  modes  would  he  accentuated  in  the  region  south  of 
the  Polar  Front,  because  all  modes  are  of  the  RSR  class. 

Figure  31  gives  the  corresponding  prediction  for  an  18-meter  SUS  on  this 
radial.  It  is  interesting  that  the  average  signal  level  increases  with  range,  the 
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Figure  21.  Radial  SW:  Parabolic  equation  propagation  loss  prediction  for  18-meter  SUS,  125  Hz  frequency. 

level  at  2000  km  being  some  10  to  15  dB  greater  than  that  at  500  km.  In 
addition,  we  see  that  the  convergenee  /one  structure  abruptly  disappears  at 
the  Polar  Front. 

As  a final  example  of  changes  in  the  propagation  characteristics  induced  by 
different  oceanographic  conditions,  the  predictions  at  63  11/  for  an  18-meter 
and  a 244-meter  SIJS  on  radial  NW  are  shown  in  Fig.  22  and  23,  respectively. 

These  predictions  were  initialized  by  the  ray  starter  and  run  under  the  same 
conditions  as  discussed  for  Fig.  10.  It  is  interesting  that  the  rather  subtle  change 
in  profile  shape  at  about  1000  km  (Fig.  2),  caused  by  a near-surface  trend  to  warmer 
temperatures,  could  have  such  a dramatic  effect  on  the  acoustics:  For  the  18-meter 
SUS  we  find  a sudden  increase  in  the  slope  of  the  convergence  zone  peaks  at  1000 
km,  which  indicates  a substantial  increase  in  t he  bottom  loss  suffered  by  the  RSR 
modes  at  this  range.  For  the  244-meter  SUS,  this  oceanographic  feature  is  marked 
by  the  sudden  appearance  of  convergence  zones.  I fleets  similar  to  these  have 
been  observed  during  a previous  propagation  trial  in  the  same  area  (Ref.  10). 
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Figure  22.  Radial  NW:  Parabolic  equation  propagation  loss  prediction  for  18-meter  SUS,  63  Hz  frequency. 


Figure  23.  Radial  NW.  Parabolic  equation  propagation  loss  predication  for  244-meter  SUS,  63  Hz  frequency. 
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CONCLUSIONS 


Ifc 


We  have  used  a parabolic  equation  model  to  study  acoustic  propagation  over  a 
wide  variety  of  typical  ocean  environments.  Bathymetric  conditions  encountered 
have  included  continental  slopes,  ocean  basins,  underwater  ridges,  and  seamounts. 

In  addition,  the  oceanographic  conditions  have  varied  from  those  typical  of  sub- 
tropical waters  in  the  north  to  antarctic  waters  south  of  the  Polar  front.  These 
conditions  have  presented  a severe  test  of  the  capabilities  of  current  acoustic  model- 
ing techniques. 

In  spile  of  difficulties  both  in  initializing  and  in  running  the  parabolic  equation 
model  on  the  continental  slope,  the  general  level  of  agreement  between  measure- 
ment and  prediction  is  encouraging.  On  a detailed  scale,  many  aspects  of  the  pre- 
diction are  unsatisfactory,  but  there  are  several  avenues  worth  exploring  which  may 
improve  the  agreement. 

Key  conclusions  resulting  from  this  work  are: 

1.  in  environments  where  the  water  depth  may  change  substantially,  such  as  near 
continental  slopes  or  underwater  features,  and  in  regions  where  the  water  depth  is 
marginal  for  RSR  propagation,  considerable  care  needs  to  be  exercised  in  the  choice 
of  vertical  sample  size,  Ar,  to  ensure  that  the  highest  order  modes  needed  are  sam- 
pled properly  at  all  ranges.  I o proceed  efficiently  anil  to  minimize  computer  run- 
ning costs,  the  model  should  automatically  adjust  A.  to  suit  the  local  conditions. 

2.  Small  changes  in  sound-speed  profile  may  induce  operationally  significant 
changes  in  acoustic  propagation  characteristics. 

3.  \coiistic  shadowing  behind  substantial  underwater  features  is  strongly  depen- 
dent on  the  nature  of  the  sediment  cover. 

4.  I he  knowledge  of  gross  sound-speed  characteristics  and  sea  floor  properties 
alone  may  not  be  sufficient  to  model  acoustic  propagation  characteristics  in  all 
ocean  regions.  Tor  example,  there  may  be  energy  loss  as  a result  of  mode  coupling 
either  from  surface  roughness  where  there  is  no  deep  sound  channel,  or  from  in- 
homogeneities  within  the  water  column  in  the  boundry  region  between  different 
water  mass  types. 
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Appendix 


APPLICABILITY  OF  SMALL-ANGLE  APPROXIMATION 


Here  consider  the  “small-angle”  approximation  used  to  derive  the  parabolic  equation 
(Eq.  6).  This  is  usually  expressed  by  the  requirement  that 


To  demonstrate  that  this  condition  is  valid  in  the  deep  sound  channel,  we  consider  the 
expression  lor  the  propagation  of  a normal  mode  under  the  parabolic  equation  approxi- 
mation. Comparing  Eq.  5 and  10,  we  obtain  the  following  expression  for  U: 


U(r,z)^UJz ) expf  K^r/Uk^) 


(A- 1 ) 


Differentiating  (A-l ) with  respect  to  r twice  we  have 


W 

dr 


(A-2) 


Substituting  this  last  equation  into  Eq.  8,  the  condition  to  be  satisfied  is: 


(A-3) 


where  we  have  used  the  well-known  (Ref.  1 3)  definition  of  equivalent  ray  angle,  0(j , to 
simplify  the  expression. 

Maximum  values  of  dQ  usually  significant  in  a deep  ocean  sound  channel  are  typically 
0.3  radian.  Thus,  for  this  type  of  propagation  the  small-angle  approximation  is  easily 
satisfied. 
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